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Abstract 

The invariant mass spectrum, forward-backward asymmetry, and lepton polarizations of the 
exclusive processes B —>■ K(K*)£ + £~ , I — fi, r are analyzed under supersymmetric context. 
Special attention is paid to the effects of neutral Higgs bosons (NHBs). Our analysis shows that 
the branching ratio of the process B — > K^i + can be quite largely modified by the effects of 
neutral Higgs bosons and the forward-backward asymmetry would not vanish. For the process 
B — > K*fi + fi~, the lepton transverse polarization is quite sensitive to the effects of NHBs, 
while the invariant mass spectrum, forward-backward asymmetry, and lepton longitudinal 
polarization are not. For both B — > Kt + t~ and B — > K*t + t~ , the effects of NHBs are quite 
significant. The partial decay widths of these processes are also analyzed, and our analysis 
manifest that even taking into account the theoretical uncertainties in calculating weak form 
factors, the effects of NHBs could make SUSY shown up. 
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1 Introduction 



The inclusive rare processes b X s £ + £ 
literatures©, |, 1 |§ g § 0, |, g 0, 0, 



= e, (j,, t have been intensively studied in 
As one of flavor changing neutral current pro- 
cesses, it is sensitive to fine structure of the standard model and to the possible new physics 
as well, and is expected to shed light on the existence of new physics before the possible new 
particles are produced at colliders. 

It is well known that invariant mass spectrum, forward-backward asymmetries(FBAs), and 
lepton polarizations are important observables to probe new physics, while the first two observ- 
ables are mostly analyzed. About lepton polarizations, it is known that due to the smallness of 
the mass of it, therefore electron polarizations are very difficult to be measured experimentally. 
So only the lepton polarizations of muon and tau are considered in literatures JTOJ, [L^, [L3|, [L4| . 
The longitudinal polarization of tau in B —>■ X s t + t~ has been calculated in standard model 
(SM) and several new physics scenarios |TIJ. For B — > X s l + l~ (/ = /x, r), the polarizations of 
lepton in SM are analyzed in Jl2| and it is pointed out that for the /i channel, the only signifi- 
cant component is the longitudinal polarization(P^), while all three components are sizable in 
the t channel. The analysis has been extended to supersymmetric models (SUSY) and a CP 
softly broken two Higgs doublet model in refs. |L3| and |14| respectively. The reference f| 
also gives a general model-independent analysis of the lepton polarization asymmetries in the 
process B — > X s t + t~ and it is found that the contribution from C 'lrlr + C 'lrrl is much larger 
than other scalar-type interactions. 

Compared with the inclusive processes B — > X s £ + £~,£ = e, /i, r, the theoretical study of 
the exclusive processes B — > K(K*)£ + £~ is relatively hard. For inclusive semileptonic decays 
of B, the decay rates can be calculated in heavy quark effective theory (HQET) JL5[]. However, 
for exclusive semileptonic decays of B, to make theoretical predictions, additional knowledge of 
decay form factors is needed, which is related with the calculation of hadronic transition matrix 
elements. Hadronic transition matrix elements depend on the non-perturbative properties of 
QCD, and can only be reliably calculated by using a nonperturbative method. The form factors 
for B decay into K^*' have been computed with different methods such as quark models ||16|| , 
SVZ QCD sum rules 0, light cone sum rules (LCSRs) |[J, [lj, [20], ^ 0. Compared to 



the lattice approach which mainly deal with the form factors at small recoil, the QCD sum 
rules on the light-cone can complementarily provide information of the form factors at smaller 
values of s. And they are consistent with perturbative QCD and the heavy quark limit. In 
this work, we will use the weak decay form factors calculated by using the technique of the 
light cone QCD sum rules and given in ||23|| . 

A upper limit on the branching ratio of B° — > K°*{i + [i~ has been recently given by CLEO 



0: 



BR(B° -> TTV+zO < 4.0 x 10" 



and they will be precisely measured at B factories, these exclusive processes are quite worthy 
of intensive study and have attracted many attentions [T23J, ^ |2JJ ^7], ^ |3UJ, [51], R21 [33], 
34], [35|, |36|]. In reference 



3|, by using improved theoretical calculations of the decay form 
factors in the light cone QCD sum rule approach, dilepton invariant mass spectra and the 
FBAs of these exclusive decays are analyzed in the standard model and a number of popular 
variants of the supersymmetric models. However, as the author claimed, the effects of neutral 
Higgs exchanges are neglected. For exclusive processes, as pointed out in |36| , the polarization 
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asymmetries of \i and r for B — > K*fi + n and B — > K*t + t are also accessible at the 



B-factories under construction. In reference [34|, the lepton polarizations and CP violating 



effects in 5 — > K*t + t~ are analyzed in SM and two Higgs doublet models. 

As pointed in refs. || £|, in two-Higgs-doublet models and SUSY models, neutral Higgs 
boson could contribute largely to the inclusive processes B — > X s £ + £~,£ — fi, r and greatly 
modify the branching ratio and FBA in the large tan/3 case. The effects of neutral Higgs in 
the 2HDM to polarizations of r in B —>■ Kt + t~ are analyzed in |05| |, and it was found that 
polarizations of the charged final lepton are very sensitive to the tan/3. 

In this paper, we will investigate the exclusive decay B — > K(K*)£ + £~ , £ = /i, r in SUSY 
models. We shall evaluate branching ratios and FBAs with emphasis on the effects of neutral 
Higgs and analyze lepton polarizations in MSSM. According to the analysis of p9| , different 
sources of the vector current could manifest themselves in different regions of phase space, 
for the very low s the photonic penguin dominates, while the Z penguin and W box becomes 
important towards high s. In order to search the regions of s where neutral Higgs bosons 
could contribute large, we analyze the partial decay widths of these two processes. Beside that 
they are accessible to B factories, our motivation also bases on the fact that to the inclusive 
processes B — > X s £ + £~~,£ = /i, r, neutral Higgs could make quite a large contributions at 
certain large tan/3 regions of parameter space in SUSY models, since part of supersymmetric 
contributions is proportional to tan 3 /3 [|] . Such regions considerably exist in SUGRA and M- 
theory inspired models |37j . We also analyze the effects of neutral Higgs to the position of the 
zero value of the FBA. Our results show that the branching ratio of the process B — > Kfi + fi~ 
can be quite largely modified by the effects of neutral Higgs bosons and the FBA would not 
vanish. Because the FBA for B — > K £ + £~ (£ = fi, r) vanishes if the contributions of NHBs 
vanish. The contributions of NHBs can be large enough to be observed only in SUSY and/or 
2HDM with large tan/3, and a non-zero FBA for B — > K £ + £~ would signal the existance of new 
physics. For the process B — > K*fi + fi~, the lepton transverse polarization is quite sensitive 
to the effects of NHBs, while the invariant mass spectrum, FBA , and lepton longitudinal 
polarization are not. For both B — > Kt + t~ and B — > K*t + t~ , the effects of NHBs are quite 
significant. Our analysis manifest that even taking into account the theoretical uncertainties in 
calculating weak form factors, the effects of NHBs could show SUSY up. In brief, our analysis 
manifest that effects of NHBs is quite remarkable in some regions of parameter space of SUSY, 
even for the process B — > Kji + ^ . 

The paper is organized as follows. In the section 2, the effective Hamiltonian is presented 
and the form factors given by using light cone sum rule method are briefly discussed. Basic 
formula of observables are introduced in section 3. Section 4 is devoted to the numerical 
analysis. In section 5 we make discussions and conclusions. 



2 Effective Hamiltonian and Form Factors 

By integrating out the degrees of heavy freedom from the full theory, MSSM, at electroweak(EW) 
scale, we can get the effective Hamiltonian describing the rare semileptonic decay b — > s£ + £~: 

Art 10 10 

H ee = ~WE«M + E^)Q,M) , (2.1) 
v2 i=1 i=1 
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where the first ten operators and Wilson coefficients at EW scale can be found in @, |3"8]f|, 
and last ten operators and Wilson coefficients which represent the contributions of neutral 
Higgs can be found in [|J. 

With the renormalization group equations to resum the QCD corrections, Wilson coeffi- 
cients at energy scale fi = m?, are evaluated. Theoretical uncertainties related to renormalization- 
scale can be substantially reduced when the next-leading-logarithm corrections are included 



The above Hamiltonian leads to the following free quark decay amplitude 
G pot 



M{h ^ si +i-) = ~^V*V tb {C 9 cS [s lfl Lb] 



t-ft + C l0 [s lfl Lb} t^ld 



-2m h C 7 eS 



si(Juv—Rb 
s 



hn +c Q1 [sRb] a +c Q2 [sRb] 



where is defined as [f40l [4l 



3tt 



a" 



Vi=0(ls),...,i/>(6s) 



T(V t ^i + i-)m Vt 
'm v 2 - sm B 2 - im v V Vi 



(2.2) 



(2.3) 



where s = s/m 2 b ,s=q 2 , C(fi) = (3 d + C 2 + 3 C 3 + C 4 + 3 C 5 + C 6 ), and 



Y(fi,s) = g(m c ,s)C(ii) 



-~g(l, s) (4C 3 + 4 C 4 + 3 C 5 + C 6 ) - ~<?(0, s) (C 3 + 3 C 4 ) 



--(3C 3 + C 4 + 3C 5 + C 6 ) 
y 



(2.4) 



where the function g(rh c , s) comes from one loop contributions of four-quark operators and is 
defined by 



Jl-^(2 + ^ 



4, 2 8 16 z 2 
--Inz 2 H 1 , 

9 27 9 s 1(2 + *£)arctan( 



ln( — v , l - it; 



1 



, 4z 2 < S 



4z 2 > s 



(2.5) 



The last terms in fl2.3|) are nonperturbative effects from (cc) resonance contributions, while the 
phenomenological factors Ki can be fixed from the processes [23] B — > K^Vi — > i^^f+f - and 
as given m the Table, g. 

Exclusive decays i? — > (if, K*)£ + £~ are described in terms of matrix elements of the quark 
operators in Eq. (p.2|) over meson states, which can be parametrized in terms of form factors. 



1 In our previous papers, e.g., || m, we follow the convention of ref. |IJ] for the indices of operators as well 
as Wilson coefficients. In this paper, we use more popular conventions (see, e.g., [|9)). That is, Os — > Og and 
Og -> Oio- 
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For the process B — > K£ + £ , the no n- vanishing matrix elements are (q = p B — p) 

(K(p)\s llx b\B(p B )) = f + (s) (p B + p), - -* * q, + -* * f (s) q„ (2.6) 



and 

(K(p)\sa^(l + ^ 5 )b\B(p B )) = (K(p)\sa^b\B(p B )) 

= *{(PB+p),s-q,(m 2 B -m 2 K )} J^^ - 
While for B — > K*£ + £~, related transition matrix elements are 

Ms) 



(2.7) 



(K*(p)\(V-A) f ,\B(p B )) = -iel(m B + m K *)A 1 (s)+i(p B +p) lx (e*p B ) 



m B + m K * 



2itlk* 



(A 3 (s) - A (s)) + e» upa e*»p B p° 



(2.8) 



and 



(X*|^^(l + 75 )6|S(p B )) = ie^e^PsP^T^s) 

+ T 2(s) {e* (m 2 B - m 2 K ,) - (ep B ) (pb + p)^} 



+ T 3 (s)(e*p B ) - 



m B - m" K * 



(Pb+p)A (2.9) 



Where e M is polarization vector of the vector meson K*. By means of the equation of 
motion, one obtains several relations between form factors 

m B + m K * m B - m K * 
M s ) = — o M s ) o M s )> 



Ivtik* 
Ao(0) = MO), 
{K*\d^\B) = 2m K *(e*p B )A (s), 

71(0) =T 2 (0). 



Ivtik* 



(2.10) 



K 


J/tf 




K 

K* 


2.70 
1.65 


3.51 
2.36 



Table 1: Fudge factors in B 
form factors. 



tfWj/W,*' _> K(*k + r decays calculated using the LCSR 
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All signs are denned in such a way as to render the form factors real and positive. The physical 
range in s extends from s min = 4m| to s max = (1 — rhx^K*) 2 - 

The calculation of the form factors given above is a real task, and one has to rely on certain 
approximate methods. We use the results calculated by using technique of LCSRs and given 



in 1 23]. And form factors can be parametrized as 

F(s) = F(0) exp( Cl s + c 2 s 2 + c 3 s 3 ). (2.11) 

The parameterization formula works within 1% accuracy for s < 15GeV 2 and can avoid the 
spurious singularities at s = m 2 B . Related parameters is given in the Table. 4 of |23| 



3 Formula of Observables 

In this section we provide formula of experimental observables, which include dilepton invariant 
mass spectrum, FBA, and lepton polarizations. 

^From eqs. (2.2 - 2.8), it is straightforward to obtain the matrix element of B —>■ K(K*)l + l~ 
as follows. 

M = - ^ V*V tb m B [T; (l>f £) + T 2 (tr 75 1) + S («)] , (3.1) 
where for B -> K£+£~, 

Tl = A!{s)p„ (3.2) 

= C'(s)p, + D'(s)q,, (3.3) 

S = 5i(a) (3.4) 

and for B -> K*£+£~, 

T; = A(s) e, pa pe* p p B p^ - iB(s) e\ + iC{s) (e* ■ p B )p, , (3.5) 

= E(s) e^e* p p a B p^ - iF(s) e\ + iG{s) (e* • p B )p„ + iH{s) (e* • p B )q^ , (3.6) 

S = i2m K *(e* ■ pb)S 2 (s) (3.7) 

with p = p B + pk,k*- Note that, using the equation of motion for lepton fields, the terms in 
in T^ 1 vanish. 

The auxiliary functions above are defined as 

A'(s) = C 9 cS (s) U(s) + ^p^C 7 cS f T (s) , (3.8) 

1 + m K 

C'(s) = C 10 f + (s), (3.9) 

DVs)-c wf ^)-^g^c Q M^ (3.10) 
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A{§) 
B{s) 
C{s) 

m 

F{s) 
G(s) 



1 — m 



K 



{m b - rh s ) 
2 



CqMS) 



„ Cffr 



1 + rh K 
(1 + rh K * 

1 

1 — m 2 K , 



C 9 ett {s)V{s) + —^C 7 ctt T 1 {s) 

s 



Cg^^Atis) + ^(1 - m K *)C 7 eti T 2 {s 



2m 



cSr 



1 - m K *)C 9 cS {s)A 2 {§) + 2m b C 7 ctt [ T 3 (s) + 



cff 



1 — m 



K* 



T 2 (s) 



T ——C w V(s) , 

(l+m x *)CioAi(s) 
1 



1 + 



C 



= ^ [(1 + m K ,)A 1 (s) - (1 - m^)A 2 (s) - 2m K *A {s)\ 



m K * 



&(§) 



1 



-4)(s)Cq2 



(m 6 + m s ) 



A (s)C Q1 



where 



fo(§) 



1 



.[$f_(S) + (l-m 2 K )f + {s)] 



1 — 

and to get the auxiliary functions given above, we have used equations of motion 

q^ii'iiM = (mi - m 2 )4> 1 ip 2 , 
<? M (^17^75^2) = -{m\ + m 2 )i) 1 ^ip 2 . 

The contributions of NHBs have been incorporated in the terms of S\(s), D'{s), H{s) 
S 2 {s). It is remarkable that the contributions of NHBs in D'(s) and H(s) are proportional to 
the inverse mass of the lepton, and for the case I = /x, the effects of NHBs can be manifested 
through these terms. 

A phenomenological effective Hamiltonian is recently given in [30|. If ignoring tensor type 
interactions in the phenomenological Hamiltonian (it is shown that physical observables are 
not sensitive to the presence of tensor type interactions ||), it is easy to verify that the matrix 
element of B — > K^*H + £~ can always be expressed as the form of the equation ( |3.1|) with the 
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auxiliary functions defined as 



where 



W c 9 + Wc 10 



A\s) = w Cl U(s) - * + Th ; 10 /r(g), 
C\s)=w C2 f + (s), 

D'(s) = w C2 Us) - 2^7^)^/0(3) , 

777,^ 

Ms) = T— ^ w Cl \/(s) - ^-(w Cg +w Cl0 )T 1 (s) , 

1 + itlk* s 



B(s) = -(1 + m K * 



w^A^s) + 4(1 - m K *)(w C9 + w clQ )T 2 (s) 
s 



1 



c ( s ) = ~1 K 1 - rh K *)w C3 (s)A 2 (s) 



/ 1 _ ^2 1 

+ (w C9 - iu ci0 ) (1 + m^-.)T 3 (s) H z J ^T 2 (s) 



F(s) = -(1 +m^*)w C4 Ai(s) , 

G(s) = — w CA A 2 (s) , 

1 + rriK* 

2rriK* , ,,,, m^-* 

= —w C4 {A 3 (s) - A (s)) + — — — — -w cg A (s) , 

s m e (m b + m s ) 



w ci = \{Cll + C LR + C RL + C RR ) , 
u! C2 = -^(—Cll + Clr — Crl + Crr) , 



w C3 — -(—Cll — Clr + Crl + C RR ) , 



w C4 


= \(c LL - 


Clr, — Crl + Crr, 




(3.38) 






+ CrlLR + ClRRL 


+ Crlrl) , 


(3.39) 




= -^{Clrlr 


+ Crllr — Clrrl 


— Crlrl) , 


(3.40) 


w C7 


= -^{Clrlr 


— Crllr + Clrrl 


~ Crlrl) , 


(3.41) 




= -^{Clrlr 


— Crllr — Clrrl 


+ Crlrl) , 


(3.42) 


Wc 9 


= rribCBR , 






(3.43) 


w c 10 


= rn s C SL ■ 






(3.44) 



In the above equations Cll,Clr, etc are defined in ref. ||. Therefore our formula given 
below can also be used to make model independent phenomenological analysis, if using Eqs. 
((ggp-dOg)) instead of Eqs. ( (|378|) - (|3TT9|) ) . 

Keeping the lepton mass, we find the double differential decay widths T K and T K * for the 
decays B — > K£ + £~ and B — > K*£ + £~, respectively, as 

a i Lx F a m B , T/ * T/ ,2 

V ts V tb\ 



dsd-u 



2 n 7T 5 



x {{\A'\ 2 + \Cf )(A - u 2 ) + \S x \ 2 {s - Am 2 ) + J Re(5 1 A'*)4m,u 

+ \C'\ 2 4m e 2 (2 + 2m 2 K - s) + Re(C'D'*)8m e 2 (l - m 2 K ) + \D'\ 2 4m/s} , (3.45) 



d 2 r** G 2 F a 2 m Bur ^ r |2 



dsdu 



x • 



+ 



2 n 7T 5 

( \A\ 2 



4 
1 

Am 2 K , 



s(A + u 2 ) + 4m 2 A) + i-L (s(A + u 2 ) - 4m 2 A) + \S 2 \ 2 (s - 4m 2 )A 



|B| 2 (A - u 2 + 8m^,(s + 2m])) + |F| 2 (A - u 2 + 8m^,(s - 4m 



- 2su \Re(BE*) + Re(AF*)l + ^^[R e (S 2 B*)(s + m 2 K , - 1) + Re(<S 2 C*)A] 



+ 



A 



Am\, 



\C\\\ - u 2 ) + |Gf (A - u 2 + 4m 2 (2 + 2m^ - s 



2mfr» 



Re(5C*)(l - m|, - s)(A - m 2 ) 
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+Re(FG*) ((1 - m 2 K , - s){\ - u 2 ) + Am]\) 



ml, 



-A 



Re(FH*) - Re(GH*)(l - m 2 K .) + \H\ 2 -^s\ . 



ml, 



in 



K* 



Here the kinematic variables (s, ii) are denned as 

s = q 2 = (p + +p_) 2 



which are bounded as 



U = {PB-P-Y ~ (PB-P+)' 

(2m,) 2 < a < (l-m K , K *) 2 
—u(s) < u < u(s) , 



with in? = me/mB and 



uls) 



'A(l-4 



mp 



A = 1 + m% K , + s 2 -2s- 2m 2 KK ,(l + s) 



V 



1- 



4m 2 



(3.46) 

(3.47) 
(3.48) 

(3.49) 
(3.50) 

(3.51) 
(3.52) 

(3.53) 



Note that the variable u corresponds to 8, the angle between the momentum of the 5-meson 
and the positively charged lepton £ + in the dilepton CMS frame, through the relation u = 
-u(s)cos8 |4g. 

Integrating over u in the kinematic region given in Eq. ( 3.50Q we get the formula of dilepton 
invariant mass spectra (IMS) 



dr- _GWm% KVtbfmD K 



ds 



D 



K 



2 10 7T 5 



(\A'\ 2 +\C'\ 2 )(X 



(3.54) 



U(S 



|Si| 2 (s-4m 2 



+ |C"| 2 4m/(2 + 2m 2 K - s) + Re(C'D' )8m/(l - m 2 K ) + \D'\ 2 4mis , 
dV^_ GWrrA > a , 

ds ~ 2 10 7T 5 ' tS [ ] 

K " = + 1% + \E\ 2 s^f- + \S 2 \ 2 (S - 4m 2 )A 

o S o 



(3.55) 
(3.56) 



D 



4m 2 K , 



\ "(*) 2 , o^,2 . o~2^ . |c,|2/\ "(^) 2 , i.;,2, 



|5| 2 (A 



+ 8m^.(s + 2m|)) + \F\\\ 



+ 8m 2 K *(s-4mj)) 
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2m\, 



u(s 



Re(BC*)(\-^)(l-m 2 K *-s) 



+ Re(PG*)((A 



-)(1 -fflp - s) +4m^A) 



7YI j^* 



Re(FH*) - Re(GH*)(l - m\,\ 



+ -rJ-sX\H\ 



in 



(3.57) 



A"* 



Both distributions agree with the ones obtained in |23], [Jo), if Cq 1; % are set to zero. 
The differential FBA is defined as 



A FB (s) 



«(*) , dr /-o dr 

dz—— + / du—— 

o ascra J-u(s) dsdu 

W) To ' dT - 

dz-—: — h / du—— 

o asaw ascra 



For B —>■ K£ + £ decays it reads as follows 

-D ,v = -2m i u(s)Re(S 1 A l 



A AK 



dB 

For 5 — > K*£ + £~ decays it reads as follows 

a ak* ( 

D K * = u(s)l s [Re(BE*) + Re(AF*)] 



(3.58) 



+ 



mi 
m K * 



Re(S 2 P*)(l - s - m\,) - Re(S 2 C*)\ 



(3.59) 



We can read from ( [3.58D , the FBA of the process B — > K£ + £~ does not vanish when the 
contributions of NHB are taken into account. With it, our analysis below also show that the 
contributions of NHBs can even be accessible in B factories. 

The lepton polarization can be defined as follows 



dT(n) 
ds 



iir, 

2 W 



1 + (Pl e L + P N e N + P T &r) ■ n 



(3.60) 



where the subscript "0" corresponds to the unpolarized width, and Pl, Pt, and P/v, correspond 
to the longitudinal, transverse and normal components of the polarization vector, respectively. 
For the process B — > K£~£ + , the Pjf , P? , and P^ , are derived respectively as 



P?D K = ^D^XRe(A'C'*) - 3m, {1 - nr h - )l{c(C""S l ) - 3m, *R<>( D' \<?, : 



(3.6i; 
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p k d k = ^Vsu(s) f _ lm{A , s ^ + 2 fi leSm {C'D' 



(3.62) 



+sRe(A'P/ 



(1 - m^)Re(A'C 
(s — 4m|) 



Re(C'5r; 



(3.63) 



D K is defined in Eq. (g3g) . For the process £ -> K*£~£ + , the Pf *, Pf *, and Pf *, are 
derived respectively as 

P**P>** = p{^ Re ( A P*) + (A + l 2 2 ^ } Re(PP*) 



A(l - m^, - s) 



A 2 



-Re(PG* + CP*) + 5— Re(CG* 

3m^* 3?7%* 



+ 



2rh e \ 



m K * 



Re(P<S 2 *) - sRe(P<S 2 ) - (1 - m 2 K *)Re(GS;) 



(3.64) 



r N u 



iry su(s) J 



AlflK 1 T^R"* 



Im(PG*)(l + 3m 2 Y » - s) 



+ Im(PPP)(l - - s) - \m{GH*)\ 



2m K *m e [Im(PP*) + Im(AF*)] 



PfD K * 



- (1 - m 2 K , - s)lm(BSZ) + AIm(C<S 2 *) 



(3.65) 



7T 



^LsRe(AB*) 



+ 



'1 — m|-« — s 



m 



K 



+ 



A 



rn 



K 

s — 4m 2 , 



+ 

m K *mi 
D K * is defined by Eq. fFBl) . 



•Re(PP*) + (1 - m^)Re(PG*) + sRe(BH* 
Re(CP*) - (1 - m|»)Re(CG*) - sRe(CPP) 
1 - m 2 ^ - s)Re(PS 2 ) - ARe(GS 2 )l 1 



(3.66) 
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nib 


^.O Vorti V 


m c 


l.^t VjtJ V 


IIO s 


n 9 GeV 

w . — v.' v V 




U.ll ^Jt3 V 


Tfl 


1 78 GeV 

1. 1 U VJC V 


iv±b 


c; os Cp\T 

O.ZjO vjc V 


1V±k 


n AQ PpV 


IVlK* 


fl 8Q PpV 

u.oy \jh v 


A/T , , , ( A/T , ,\ 


q in/''! fiCA PoV 


r R 


4 22 x 10" 13 GeV 




8.70(27.70) x 10~ 5 GeV 


r(j/v> -> 


5.26 x 10~ 6 GeV 




2.14 x 10~ 6 GeV 




1.17 x 10~ 5 GeV" 2 


a" 1 


129 


l^ 6 | 


0.0385 



Table 2: Values of the input parameters used in our numerical analysis. 





c 2 


c 3 


c 4 


c 5 


c e 


cf 


c 9 


Cio 


C 


-0.248 


+ 1.107 


+0.011 


-0.026 


+0.007 


-0.031 


-0.313 


+4.344 


-4.669 


+0.362 



Table 3: Wilson coefficients of the SM used in the numerical analysis. 

4 numerical analysis 

Parameters used in our analysis are list in Table |[ Considering that the branching ratios of 
B — > K £ + i~ and B — > K* £ + £~ are not very sensitive to the mass of m&, we neglect the 
difference between the pole mass and running mass of b quark. 

The Wilson coefficients in the SM used in the numerical analysis is given in the Table |||. 
Cy-^ is defined as 

C e 7 ff = C 7 - C 5 /3 - C 6 . (4.67) 



Cq 1 2 come from exchanging NHBs and are proportional to tan 3 /? in some regions of the 
parameter space in SUSY models. According to the analysis in ||, |37j , the necessary conditions 
for the large contributions of NHBs include: (i) the ratio of vacuum expectation value, tan/?, 
should be large, (ii) the mass values of the lighter chargino and the lighter stop should not be 
too large (say less than 120 GeV), (ra)mass splitting of charginos and stops should be large, 
which also indicate large mixing between stop sector and chargino sector. As the conditions 
are satisfied, the process B — > A s 7 will impose a constraint on C7. It is well known that 
this process puts a very stringent constraint on the possible new physics and that SUSY can 
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SUSY models 


R 7 


Rg 


-Rio 




Cq 2 


SUSY I 


-1.2 


1.1 


0.8 


0.0 


0.0 


SUSY II 


-1.2 


1.1 


0.8 


6.5(16.5) 


-6.5(-16.5) 


SUSY III 


1.2 


1.1 


0.8 


1.2(4.5) 


-1.2(-4.5) 



Table 4: Wilson coefficients of the SUSY used in our numerical analysis. Ri means Ci/Cf M . 
SUSY I corresponds to the regions where SUSY can destructively contribute and can change 
the sign of C-j, but the contributions of NHBs are neglected. SUSY II corresponds to the regions 
where tan/3 is large and the masses of superpartners are relatively small. SUSY III corresponds 
to the regions where tan/3 is large but the masses of superpartners are relatively large. In the 
last two cases the effects of NHBs are taken into account. The contributions of NHBs are 
settled to be different for both the case £ = fi and £ = r, since Cq x 2 are proportional to the 
mass of lepton. The values in bracket are for the case £ = r . 

2 

contribute destructively when the signature of the Higgs mass term /x is minus. There exist 
considerable regions of SUSY parameter space in which NHBs can largely contribute to the 
process b — > s£ + £~ while the constraint of b — > 57 is respected (i.e., the signature of the Wilson 
coefficient C7 is changed from positive to negative). When the masses of SUSY particle are 
relatively heavy (say, 450 Gev), there are still significant regions in the parameter space of 
SUSY models in which NHBs could contribute largely. However, at these cases C7 does not 
change its sign, because contributions of charged Higgs and charginos cancel with each other. 
We will see it is hopeful to distinguish these two kinds of regions of SUSY parameter space 
through observing B — > K^*H + £~ . 

As pointed out in || |J, the contribution of NHBs is proportional to the lepton mass, 
therefore for £ = e, contributions of NHBs can be safely neglected. While for cases £ — fi and 
£ = t, the contributions of NHBs can be considerably large. To investigate the effects of NHBs 
in SUSY models, we take typical values of 67,9,10 an d Cq 12 as given in Table 6. The SUSY 
model without considering the effects of NHBs (SUSY I in Table 4) is given as a reference 
frame so that could the effects of NHBs be shown in high relief. 

Numerical results are shown in Figs. 1-4. In Fig. |l|(a), the IMS of B — > Kfi + fi~ is depicted. 
We see that at the high s regions, NHBs greatly modify the spectrum. While at the low s 
region, the effects of NHBs become weak. In Fig. |l|(b), the FBA of the B — > K[i + fi~ is 
presented. Fig. |l](b) shows that the average FBA in B —>■ K{i + fi~ is 0.02. To measure an 
asymmetry A of a decay with the branching ratio Br at the na level, the required number 
of events is iV = n 2 /(BrA 2 ). For B — > Kfi + fi~, the average FBA is 0.02 or so, the required 
number of events is 10 12 or so. Therefore it is hard to observe the derivation of FBA from 
the SM. In Fig. |l](c) and Fig. |l|(d), the longitudinal and transverse polarizations are given. 
The effect of NHBs to the longitudinal polarization is weak but the effect to the transverse is 
remarkable. 

In Fig. 0(a) and Fig. 0(b) the IMS and FBA of B — > Kt + t~ are presented respectively. 
For SUSY II, the effects of NHBs to IMS is quite manifest, and the average FBA can reach 
0.1. For SUSY III, the average FBA can reach 0.3. Therefore, in order to observe FBA, the 
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(a) (b) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.' 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



(c) (d) 

Figure 1: The IMS (a), FBA(b), P L (c), and P T (d) of the process B -> Kfj+fj,-. The solid 
line, dashed line, dot line and dashed-dot line represent the SM, SUSY I, SUSY II, SUSY III 
respectively. Both the total (SD+LD) and the pure SD contributions are shown in order to 
compare. 
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(b) 







0.8 

















" ^^^^ 


0.6 






-0.2 




0.4 
















-0.4 




0.2 

















-0.6 




-0.2 






-0.B 




-0.4 










-0.6 




I 


-1 




-0.8 















0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 O.B 0.85 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 



(c) (d) 

Figure 2: The IMS(a), FBA(b), P L (c), and P T (d) of the process B -> Kt+t' . The line 
conventions are the same as given in the legend of Fig 1. 
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required number of events should be 10 9 or so and 10 8 , respectively, so that in B factories, 
say LHCB, these two cases are accessible. In Fig. 0(c) and Fig. ||(d), the longitudinal and 
transverse polarizations are drawn respectively. The effects of NHBs are also very obvious. 

Figs. 3 and 4 are devoted to the decay B -> K*l + l~. In Fig. |, the IMS, FBA, and 
polarizations of B — > K*fi + fi~ are given. We see that this process is not as much as sensitive 
to the effect of NHB as B — > Kfi + {i~ . However, the contribution of NHBs will increase the part 
with positive FBA and will be helpful to determine the zero point of FBA. Fig. |3](d) depicts 
the transverse polarization of the B — > K*fi + fi~ , and the effect of NHBs is quite obvious. The 
zero point of the FBA can be slightly modified as shown in Figure 0(b) due to the contributions 
of NHBs. 

In Fig. |], the IMS, FBA, longitudinal and transverse polarizations of the B — > K*t + t~ 
are depicted. The effect of NHBs does show in great relief. It is worth to note that IMS, FBA, 
and lepton polarizations for B — > K* £ + £~ in MSSM without including the contributions of 
NHBs are also significantly diffident from those in SM, while for B — > K £ + £~ they have little 
differences from those in SM. Therefore, compared to the process B — > K £ + £~, more precise 
measurements for B — > K* £ + £~ are needed in order to single out the contributions of NHBs. 

Normal polarizations for both B — ► K £ + £~ and B — > K* £ + £~ are small and can be 
neglected because the imaginary parts of Wilson coefficients are small in SUSY models without 
CP violating phases which are implicitly assumed in the paper. 

The behavior of IMS (a), FBA(b), -Pl(c), and Pr(d) shown Figs 1-4 can be understood 
with the formula given in the Section 3. With Eqs. (|3.55|) , (3.10) and (3.11), we see that the 
contributions of NHBs are contained in the terms of Si and D' . At the high s regions, it is 
these two terms which are important. This explained the behavior of IMS given in (a) of Fig. 
1 and Fig. 2. The Eq. ( |3.58| ) shows that the FBA is proportional to the mass of the lepton. 
For the case B — > Kfi + fi~, due to smallness of the mass /i, the FBA does not vanish but is 
hard to be measured. While for the case B — > Kt + t~ , the mass r is quite large and observing 
FBA is relatively easy. For SUSY II, though the numerator of FBA is comparatively large, the 
large IMS suppresses the value of FBA; for SUSY III, the numerator is relatively small, but 
the FBA do demonstrate the effects of NHBs more manifestly, as shown in Fig. 2(b) due to 
smallness of IMS. The Eqs. (3.63) and (3.64) show that for the case £ — fi, the contributions of 
NHBs to P N , P T are suppressed by the mass of \i. But for the case £ = r, the contributions of 
NHBs become quite manifest both for SUSY II and SUSY III. The term with D' in Eq. (|3.63|) 
will change its sign when there exists relatively not too small contributions of NHBs, the fact 
deduced from Eq. (3.10), that explains why the sign of Pt is changed. The difference between 
the case SUSY II and SUSY III is small, the reason is just the same as stated in the analysis 
of FBA. 

Since the terms incorporating the contributions of NHBs is proportional to A as shown in 
Eq. (|3.57 ), which approaches zero at high s regions; while at small s regions, the effects of 
NHBs are dwarfed by the other contributions. Therefore, only when Cq 1 are quite large could 
effects of NHBs be manifest, as shown in Fig. 3(a) and Fig. 4(a). According to the Eq. (|3.59 ), 
at high s regions, the effects of NHBs would be suppressed by A and 1 — s — m 2 K ,. The same 
suppression mechanism exists for P^. This suppression mechanism explains the fact that the 
processes B — > K* £ + £~ are not sensitive to the effects of NHBs. However, when there exist 
large contributions of NHBs, the sign of Pt will be changed, as indicated in both Fig. 3(d) 
and Fig. 4(d). 
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(b) 





0.1 0.2 0.3 0.4 0.5 0.6 



0.1 0.2 0.3 0.4 0.5 0.6 



(c) 



(d) 



Figure 3: The IMS(a), FBA(b), P L (c), and P T (d) of the process B -> K*fi+fi~ 
conventions are the same as given in the legend of Fig 1. 



The line 
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model 


A 


B 


c 


D 


E 


tot(SD) 


tot(SD+LD) 


SM 


LCSR 


0.353 


54.707 


0.032 


4.566 


0.076 


0.573 


59.736 


SVZ 


0.215 


22.918 


0.015 


1.593 


0.026 


0.299 


24.767 


SUSY I 


LCSR 


0.425 


54.723 


0.037 


4.576 


0.086 


0.675 


59.847 


SVZ 


0.179 


22.910 


0.011 


1.586 


0.019 


0.236 


24.704 


SUSY II 


LCSR 


0.556 


54.865 


0.131 


4.833 


0.849 


2.067 


61.233 


SVZ 


0.348 


23.009 


0.068 


1.726 


0.321 


1.002 


25.473 


SUSY III 


LCSR 


0.429 


54.727 


0.040 


4.584 


0.109 


0.717 


59.889 


SVZ 


0.181 


22.912 


0.012 


1.590 


0.028 


0.255 


24.723 



Table 5: Partial decay widths for B — ► /i . LCSR means the approach light-cone 
QCD sum rules, SVZ means the SVZ QCD sum rule |[L7|| . Character A means the region 

~ 2 ~ 2 -2 -2 -2 ~2 -2 

(s , (jfty -5) ), B ((j7ty - 5) , (my + 5) ), C ((m^ + 5) , (rh^ - 5) ), D((m^ - 5) , (m^ + 5) ) 
and E ((m^ + 5) , s max ). The unit is T B x 10~ 6 , which is 4.22 x 10" 19 GeV. 5 is selected to be 
0.2 GeV. 5 is normalized with Mb 



model 


A 


B 


c 


D 


E 


tot(SD) 


tot(SD+LD) 


SM 


LCSR 


0.930 


83.257 


0.141 


9.976 


0.258 


1.882 


94.562 


SVZ 


2.943 


111.278 


0.147 


7.504 


0.137 


3.639 


122.008 


SUSY I 


LCSR 


1.627 


83.402 


0.198 


10.085 


0.330 


2.915 


95.64 


SVZ 


4.517 


111.423 


0.183 


7.552 


0.149 


5.291 


123.825 


SUSY II 


LCSR 


1.178 


83.431 


0.234 


10.164 


0.352 


2.677 


95.360 


SVZ 


2.801 


111.292 


0.156 


7.525 


0.145 


3.522 


121.918 


SUSY III 


LCSR 


1.631 


83.407 


0.201 


10.092 


0.334 


2.938 


95.664 


SVZ 


4.518 


111.425 


0.184 


7.553 


0.150 


5.296 


123.830 



Table 6: Partial decay widths for B — ► K*jjL + n . Other conventions can be found in Table 5. 



model 




A' 


B' 


tot(SD) 


tot(SD+LD) 


SM 


LCSR 


1.884 


0.094 


0.132 


1.978 


SVZ 


0.659 


0.036 


0.054 


0.695 


SUSY I 


LCSR 


1.884 


0.086 


0.131 


1.970 


SVZ 


0.655 


0.025 


0.038 


0.680 


SUSY II 


LCSR 


2.022 


1.496 


1.674 


3.519 


SVZ 


0.726 


0.552 


0.637 


1.278 


SUSY III 


LCSR 


1.874 


0.094 


0.129 


1.968 


SVZ 


0.651 


0.026 


0.035 


0.677 



Table 7: Partial decay widths of B — > Kt + t~ . Character A' means (sq, (m^ — 5) ), B' means 
((ttv + 5) 2 , s max ). The unit is T B x 10~ 6 , which is 4.22 x 10~ 19 GeV. 
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model 




A' 


B' 


tot(SD) 


tot(SD+LD) 


SM 


LCSR 


4.045 


0.096 


0.183 


4.141 


SVZ 


3.029 


0.048 


0.102 


3.076 


SUSY I 


LCSR 


4.088 


0.173 


0.327 


4.261 


SVZ 


3.052 


0.072 


0.159 


3.124 


SUSY II 


LCSR 


4.148 


0.266 


0.460 


4.413 


SVZ 


3.054 


0.084 


0.167 


3.138 


SUSY III 


LCSR 


4.078 


0.168 


0.312 


4.246 


SVZ 


3.050 


0.071 


0.156 


3.121 



Table 8: Partial decay widths of B — » K*t + t . Other conventions can be found at Table 7. 



The partial decay widths (PDW) are listed in Tables. 5, 6, 7, and 8. We see that at the high s 
region, for the process B — > Kl + l~ , l=/z, r, the contributions of NHBs do show up, as expected. 
For B — > K*l + l~ , the effects of NHBs in the high s region is signifiacnt when \—r while they 
are small for l=/j,. It can be read out from these four table that the results are consistent with 
the Fig. |(a), 0(a), 0(a), and f|(a). In order to estimate the theoretical uncertainty brought by 
the methods calculating the weak form factors, we use the form factors calculated with LCSR 
and SVZ QCD sum rules (SVZ) method 0. For B -> K £+ (r , PDWs calculated with form 
factors obtained by SVZ method is 50% of those by LCSR approach; while for B — ► K* £ + £~, 
PDWs increase 100% or so. We see that at low s regions the theoretical uncertainty can reach 
from 100% to 200%. Another point worthy of mention is that the contribution of resonences 
domainate the integerated decay width, as had been pointed out in 



5 Conclusion 

We have calculated invariant mass spectrum, FBAs, and lepton polarizations for B — > K £ + i~ 
and B — > K* £ + £~ l=/i, r in SUSY theories. In particular, we have analyzed the effects of NHBs 
to these processes. It is shown that the effects of the NHBs to B — > Kt + t~ and B — > K*t + t~ 
in some regions of parameter space of SUSY models are considerable and remarkable. The 
reason lies in the mass of the r, which can magnify the effects of NHBs and can be see through 
from the related formula. The numerical results imply that there still exist possiblities to 
observe the effects of NHB in B — > Kfi + fi~ and B —>■ K*fi + fi~ through IMS, FBA and lepton 
polarizations of these processes. In particular, for B — Kfi + fi~ in the case of SUSY II, the 
partial width in the high s where short distance physics dominants can be enhanced by a factor 
of 12 compared to SM. Our analysis also show that the theoretical uncertainties brought in 
calculating of weak form factors are quite large. But the effects of NHBs will not be washed 
out and can stand out in some regions of the parameter space in MSSM. If only partial widths 
are measured, it is difficult to observe the effects of NHBs except for the decay B — > Kt + t~ . 
However, the conbined analysis of IBS, FBA, and lepton polarizations can provide usefull 
knowledge to look for SUSY. Finally, we would like to point out that FBA for B — > Kl + l~ 
vanishes (or, more precisely, is neglegiblly small) in SM and it does not vanish in 2HDM and 
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SUSY models with large tan/3 due to the contributions of NHBs. However, only in SUSY 
models and for l=r it is large enough to be observed in B factories in the near future. 
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